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ABSTRACT 

We report the detection of a new source of very high energy (VHE, E 7 > 100 GeV) 7-ray emission located 
close to the Galactic Plane, MAGIC J0616+225, which is spatially coincident with SNR IC 443. The observa- 
tions were carried out with the MAGIC telescope in the periods December 2005 - January 2006 and December 
2006 - January 2007. Here we present results from this source, leading to a VHE 7-ray signal with a statistical 
significance of 5.7 sigma in the 2006/7 data and a measured differential 7-ray flux consistent with a power 
law, described as dA/ 7 /(dAdfd£) = (1.0±0.2) x 10- n (£/0.4TeV)- 31±0 ' 3 cirrV'TeV -1 . We briefly discuss the 
observational technique used and the procedure implemented for the data analysis. The results are placed in 
the context of the multi wave length emission and the molecular environment found in the region of IC 443. 
Subject headings: gamma rays: observations — supernovae remnants — ISM:individual (MAGIC J 0616+225, 
IC443) 



1. INTRODUCTION 

IC 443 is an asymmetric shell-type SNR with a diame- 
ter of 45 arc minutes at a distance of about 1.5 kpc (IFesenI 
Il984t IClaussen et al.lll997l) . It is included in Green's cata- 
log (IGreenl 120041) . and it has a spectral index of 0.36, and 
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a flux density of 160 Jy at 1 GHz. It was mappe d in ra- 
dio with the VLA at 90 cm (IClaussen et all 1 1997b and a t 
20, 6 and 3.5 cm (lolbertetal.ll200U ICondon et all. fl998l). 
Moreover, IClaussen et al.l dl997b reported the presence of 
maser emission at 1720 MHz from four sources, the strongest 
of which is located at (l,b) ~ (-171.0,2.9). IC 443 is 
a prominent X-ray source, w ith data available from Rosa t 
(lAsaoka & Aschenbachll 19941). ASCA dKeohane et al.ll 19971) 
XMM dBocchino & BykovT HoOO. 20011 120031: iBvkoy et al l 



2005 1 iTroia et al . ._ 
Gaensler et all 12006). 



2006), and Chandra Joibert et alj I200H 
The EGRET has detected a 7-ray 



source above 100 M e V in the IC 443 SNR, 3 EG J0617+2238 
dEposito et al .1119961: iHartman et ail 1 19991) . Upper limits to 
the very high energy (VHE) 7-ray emission from IC 443 
were reported by the Whipple collaboration: dA/ 7 /(dAdf) < 
6 x Kr^cnrV 1 (0.11 Crab) above 500 GeV (lHolderll2005l) 
and by the CAT collaboration d iV-y /(dAdf ) < 9 x Kr 12 cirr 2 s _1 
above 250 GeV (lKhelifi||2003l) . 

Here we present observations of the SNR IC 443 with the 
Major Atmospheric Gamma Imaging Cherenkov (MAGIC) 
telescope resulting in the detection of a new source of VHE 
7-rays, named MAGIC J06 16+225. We briefly discuss the 
observational technique used and the procedure implemented 
for the data analysis, derive a VHE 7-ray spectrum, and an- 
alyze it in comparison with other observations, including the 
molecular environment found in the region of IC 443. 
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2. OBSERVATIONS 

MA GIC (see e.g., iBaixeras etlill (|2004|) : ICortina et all 
(2005) for a detailed description) is the largest single dish 
Imaging Air Cherenkov Telescope (I ACT) in operation 25 . It 
is located on the Canary Island La Palma (28.8°N, 17.8°W, 
2200 m a.s.l.). 

The SNR IC 443 was observed for a total of 10 hours in 
the period December 2005 - January 2006 (period I), with 
the telescope pointing to the SNR center. The VHE 7-ray 
sky map showed evidence for a VHE 7-ray signal with a sta- 
tistical significance of 3er (before trials were taken into ac- 
count). To test the hypothesis that this excess is due to a 
VHE 7-ray source, MAGIC J0616+225, the excess center 
was observed for a total of 37 hours in the period Decem- 
ber 2006 - January 2007 (period II). Therefore, in the analysis 
of the period II data no trial factors needed to be taken into 
account to compute the statistical significance of the source 
detection. Changes in the readout chain in spring 2006 due 
to the installation of a novel 2 GSam ples/s FADC system 
dBartko et alj|2005t iGoebel et aT]|2007l) in parallel to the ex- 
isting 300 MSamples/s FADC system made it advisable (on 
grounds of simplicity of the analysis) to select only the sec- 
ond part of the dataset, period II, for the studies presented 
here. However, within statistics the data from period I pro- 
vide compatible results. 

At La Palma, IC 443 culminates at about 6° zenith angle 
(ZA). The observations were carried out in the false-source 
tracking (wobble) mode dFomin et al][l9 94). The sky direc- 
tions (Wl, W2) to be tracked were on two opposite sides of 
the source direction, at a distance of 0.4° from the source. 
They were chosen such that in the camera the sky bright- 
ness distribution relative to Wl was similar to the one rela- 
tive to W2. For each tracking position three circular back- 
ground control regions were defined, located symmetrically 
to the source region with respect to the camera center. Dur- 
ing wobble mode data taking, 50% of the data was taken at 
Wl and 50% at W2, switching (wobbling) between the 2 di- 
rections every 20 minutes. In total, about 8 million triggers 
were recorded in period I and about 30 million triggers were 
recorded in period II. There are two bright stars in the field of 
view: rj Gem (mag(V) = 3.28 and B- V = 1.60) at a distance 
to MAGIC J0616+225 of 0.4° and ^ Gem (mag(V) = 2.88 
and B- V = 1 .64) at a distance of 1 .9° to MAGIC J0616+225. 
Both stars are rather red such that the increase in the pixel an- 
ode currents is still sustainable at nomin al pixel high voltage 
dAlbert et alj|2007bt IPaneque et alj|2004l) . The discriminator 
thresholds of the channels included in the trigger are dynami- 
cally regulated to keep th e individual pixel rates at a constant 
level ( ICortina et alj|2005l) . 

3. DATA ANALYSIS 

Data runs with anomalous trigger rates due to bad obser- 
vation conditions (three nights of data taking) were rejected 
from further analysis. In addition, data with ZA> 30° (2 h) 
were excluded. The remaining data set corresponded to an ef- 
fective observation time (including a dead-time correction) of 
29 h. 

The data analysis was carried out using the stan- 
dard MAGIC analy sis and reconstruction software 
dBretz & Wagnerl 12003b . the first step of which involves 
the FADC signal reconstructi on using a digital filter and the 
calibration of the raw data dAlbert et alj l2006dt iGaug et ail 

25 see also http://wwwmagic.mppmu.mpg.de/magic/factsheet7 
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FIG. 1 . — Sky map of 7-ray candidate events (background subtracted) in the 
direction of MAGIC J06 16+225 for an energy threshold of about 150 GeV 
in galactic coordin ates. Overlayed are 12 CO emission contou rs (cyan) from 
IDame et"aT] {2001), contours of 20 cm VLA radio data from | C ondon e t all 
11998) (green), X-ray contours from Rosat I Asaoka & Aschenbach 1994) 
(purple) and 7-ray contours from EGRET (Hartman et al. 1993) (black). The 
12 CO contours are at 7 and 14 K km/s, integrated from -20 to 20 km/s in ve- 
locity, the range that best defines the molecular cloud associated with IC 443. 
The contours of the radio emission are at 5 mjy/beam, chosen for best show- 
ing both the SNR IC 443. The X-ray contours are at 700 and 1200 counts / 
6 ■ 10~ 7 sr. The EGRET contours represent a 68% and 95% statistical proba- 
bility that a single source lies within the given contour. T he white star denotes 
the position of the pulsar CXOU J061705.3+222127 IQlbert et alj200ll). The 
black dot shows the position of the 1720 MHz OH maser IClaussen et all 
\\Wft) . The white circle shows the MAGIC PSF of cr = 0.1°. 
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FIG. 2. — Distribution of excess 7-ray candidate events (figure ffj 
as a function of the squared angular distance from the excess center of 
MAGIC J0616+225 (points), compared to the expected distribution for a 
point-like source (blue line) corresponding to the MAGIC PSF. 



I2005I) . Thereafter, shower images were cleaned by applying 
a cut of 10 photoelectrons (ph. el.) for image core pix els and 
5 ph. el. for boundary pixels, respectively (see e.g. dFeganl 
1997)). These tail cuts were scaled for the larger size of the 
outer pixels of the MAGIC camera. T he camera im ages were 
characterized by image parameters dHillasi 119851). In thi s 
analysis, t he Ra ndom Forest method (see iBock et alj d2004l) : 
Breiman (2001) for a detailed description) was app l ied fo r 
the 7/hadron separation (for a review see e.g. iFeganl 1997)) 
as well as for the energy estimation, see e.g. lAlbert et all 
(2006aifb|R I2007bl) . Above the analysis energy threshold, 
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FIG. 3. — VHE 7-ray spectrum of MAGIC J06 16+225 (statistical errors 
only). The full line shows the result of a simple power law fit to the spectral 
points taking into account correlations between the spectral points that are 
introduced by the unfolding procedure. The dashed line shows t he spectrum 
of the Crab Nebula as measured by MAGIC jAlbert et alJ2007bl) . 

the cut efficiency reaches about 50% corresponding to an 
effective collection area for 7-ray showers of about 50000 
m 2 . 

For each event the arrival direction of the primary 
7-ray candidate in sky coordinates was estimated using 
the DISP-method dFomin et all Il994t iLessard et ail 1200 U 
iDomin go-Santamari a et al.l l2005lh A conservative cut of 
SIZE> 200 ph. el. was applied to select a subset of events 
with superior angular resolution. For these events the effect 
of the bright stars in the field of view is negligible. The corre- 
sponding analysis energy threshold (the energy corresponding 
to the maximum of the differential 7-ray rate after all analysis 
cuts ) was about 150 GeV. 

Figure Q] shows the sky map of 7-ray candidates (back- 
ground subtracted) from the direction of MAGIC J06 16+225 
in galactic coordinates. It was folded with a two-dimensional 
Gaussian with a standard deviation of 0.072°. The MAGIC 
7-ray PSF (standard deviation of a two-dimensional 
Gaussian fit to the non-folded brightness profile of a 
point source) is (0.10 ± 0.01) ° for SIZE> 200 ph.el. 
(Domingo-Santamaria et al. 2005). The folding of the sky 
map serves to increase the signal-to-noise ratio by smoothing 
out statistical fluctuations. However, it correlates the values 
of neighbouring bins and somewhat degrades the spatial reso- 
lution. The sky map is overlayed wi th contours of 12 CO emis- 
sion (cyan) from Dame et al. (2001), contours of 20 cm VLA 
radio data from Condon et al. ( 1998) (g reen), X-ra y contours 
from Rosat (purple) ( Asaoka & Aschenbach 1994) and 7-ray 
contours from EGRET dHartman et al.l 1 1999) (black). Also 
XMM-Newton and Chandra X-ray observations of IC 443 ex- 
ist; they will be discussed below. The white sta r denotes the 
positi on of the pulsar CXOU J06 1 705 .3+222 1 27 dOlbert et alj 
1200 11) . Th e black dot shows the position of the 1720 MHz 
OH maser (IClaussen et alJfl997l) . The VHE 7-ray sky map 
shows a clear excess centered at (RA, DEC)=(06 h 16 m 43 s , 
+22°31'48"), MAGIC J0616+225. The statistical position 
error is 1.5', the systematic error due to background de- 
termination and poi nting uncertainty is estimated to be 1' 
dAlbertetalj|2007bl) . Within errors MAGIC J0616+225 is 
point-like. The center of gravity of the excess in the period 
II data presented here agrees with the center of gravity of the 
excess in the period I data not shown here. 

Figure [2] shows for the excess 7-ray-like events in figure 
Q] the distribution of the squared angular distance, 9 2 , with 
respect to the excess center together with the expected distri- 



bution for a point-like source. Chosing a conservative cut 
of 9 2 < 0.05 deg 2 , as appropriate for unidentified sources, the 
observed excess in the period II data set in the direction of 
MAGIC J06 16+225 has a significance of 5.7a. The period I 
data showed an excess of about 3a significance at a position 
less then 3' = 0.05° away from period II data set. Computing 
the significance of the period II data set with respect to the 
period I position yields only a negligible difference for the 
chosen large signal region 2 < 0.05 deg 2 . Chosing smaller 2 
cuts, one obtains a higher significance for the period II excess 
center. 

For the spectral analysis the excess data from a sky re- 
gion of maximum angular distance of 2 = 0.05 deg 2 around 
the excess center were integrated. Figure [3] shows the re- 
constructed VHE 7-ray spectrum (dN 7 /(dE 7 dAdt) vs. true 
E 7 ) of MAGIC J06 16+225 after correcting (unfolding) for the 
instru mental energy resolution ( Anvke ev et alJ[T99 lh iBerterol 
1989). The horizontal bars indicate the bin size in energy. A 
simple power law was fitted to the spectral points taking into 
account correlations between the spectral points that are in- 
troduced by the unfolding procedure. The result is given by 
(X 2 /n.d.f=l.l): 



- (L0±0 - 2)Xl ° (o4Tw) 



dAdtdE 



-3.1+0.3 



cnrV'TeV" 1 . 



The quoted errors are statistical. The systematic error is esti- 
mated to be 35% in the flux level determination an d 0.2 in the 
spectral index, see also lAlbert et al.l (2006b. 2007b). The inte- 
gral flux of MAGIC J0616+225 above 100 GeV is about 6.5% 
and above 300 GeV about 2.8% of the Crab Nebula flux. The 
integral flux of MAGIC J06 16+225 in the observation period 
II is compatible within errors with that of period I. Within the 
two observation periods (two months each) no flux variations 
exceeding the measurement errors were observed. All data 
analysis steps were cross-checked by a second, independent 
analysis, yielding compatible results. 

4. DISCUSSION AND CONCLUDING REMARKS 

If, due to the spatial association shown in Figure Q] we ac- 
cept that MAGIC J0616+225 is associated with SNR IC 443, 
located at a distance of ~ 1.5 kpc, then it has a luminosity 
between 100 GeV and 1 TeV of about 2.7 x 10 33 erg s _1 . This 
is roughly one order of magnitude below t he luminosity of 
the so urces HESS J1813-178 and J183 4-087 dAharqnian et al] 
120061) . both also detected by MAGIC dAlbert et al Il2006allcl) . 
In addition, the 7-ray spectrum of MAGIC J06 16+225 is 
steeper than those of J1813-178 and J1834-087. 

Figure Q] indeed shows a very interesting multi-frequency 
phenomenology. First, it can be noted that the MAGIC VHE 
7-ray source is slightly displaced with respect to the central 
position of the EGRET source 3EG J0617+2238, although 
still consistent with it, within errors. An independent anal- 
ysis of GeV photons me asured by EGRET resu lted in the 
source GeV J06 17+2237 dLamb & Macombll 19971). th at is at 
the same location o f 3E G J06 17+223 8. iGaisser et all d!998l) : 
Barin g et ail d 19991) and iKaull d2001l) extrapolate the energy 
spectrum of 3EG J0617+2238 into the VHE 7-ray range. The 
predicted fluxes are higher with a harder energy spectrum than 
MAGIC J06 16+225. This supports the view that a direct ex- 
trapolation of the EGRET source to the MAGIC energy range 
is not valid. 

The EGRET source is located in the center of the SNR, 
whereas the VHE 7-ray source is displaced to the south, 
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in direct correlation with a molecular cloud. The molecu- 
lar cloud environment surrounding IC 443 and its possible 
connection with th e EGRET 7-ray source was studied by 
iTorres et"aT] d2003l) . who also provided a review of previ- 
ous measurements regarding this SNR environment. There 
is a large amount of molecular mass (< 10 4 M©) consis- 
tent with the distance to the SNR IC 443, corresponding to 
a velocity range of -20 to 20 km/s, as shown in Figure Q] 
The highest CO intensity detected is directly superimposed 
on the central positio n of the MAGIC source. Moreover, 
IClaussen et alj d 1997b reported the presence of maser emis- 
sion from (l,b) ~ (-171.0 , 2.9), spatial l y corr elated with the 
MAGIC source (see also Hewitt et ail (I2006I) '). The maser 
emission is an indication for a shock in a high matter density 
environment. It is assumed to be due to collisionally excited 
H2 molecules heated by the shock. 

An electronic bremsstr ahlung hypothesis f or the origin of 
the EGRET source (e.g., iBvkov et alj (120001) ) is difficult to 
reconcile with the fact that the radio synchrotron, X-rays, and 
optical emission is concentrated towards the rim of the rem- 
nant, whereas the EGRET source is located in the center. On 
the other hand, the optical emission seems to fade in regions 
where CO e mission increases an d where the MAGIC source is 
located (see lLasker et al.l dl990l) ). This perhaps indicates that 
the molecular material, which absorbs the opti cal radiation, is 
at the remnant's nearest side. The report by ICorneta (1 19771) 
also argues that the molecular mass is located between us and 
th e SNR. Simi l arly, t he recent analysis of XMM observations 
by lTroja et al.l (|2006|) reached the same conclusions. The ob- 
served VHE 7-radiation may be due to 7r°-decays from in- 
teractions between cosmic rays accelerated in IC 443 and the 
dense molecular cloud. A possible distance of this cloud from 
IC 443 could explain the steepness of the VHE 7-ray spectrum 
measu red. As has been emphasized by I Aharonian & Atoyanl 
( 1996), the observed 7-rays can have a significantly different 
spectrum from that expected from the particle population at 
the source (the SNR shock). 

The positions of 3 EG J0617+2238, GeV J0617+2237, and 



MAGIC J06 16+225 are all different from that o f the pulsar 
wind nebula (PWN) CXOU J06 1705.3+222127 (lOlbert et alJ 
l200UlBocchino & Bvkovl2001l) . The PWN is now co-located 
with a high density molecular material region dSeta et ail 
1998), which in addition is excited, as measured by a high 
CO(J = 2- l)/CO(7 =1-0) ratio. If the VHE 7-ray emis- 
sion were related to the PWN, one may expect some spatial 
overlap between the PWN and the 7-ray sources. 

A complete coverage of the X-ray emission from the region 
was made with XMM dBocchino & Bvkovll2003l) . resulting 
in the detection of 12 X-ray sources with fluxes larger than 
5 x 10~ 14 erg cm" 2 s" 1 . None of these sources is spatially 
coincident with the MAGIC detection reported here. Rather, 
they are mostly located in the relatively small, 15 arcmin x 
15 arcmin region, for which the analysis of the 2MASS data 
reveals strong 2.2 /im emission indicating interaction with a 
molecular cloud. The MAGIC source, uncorrelated with X- 
ray sources, is, however, also co-spatial with a region of high 
2.2 /im emission, but farther away from the shock. 

In summary, the observation of IC 443 using the MAGIC 
Telescope has led to the discovery of a new source of VHE 
7-rays, MAGIC J0616+225, near the Galactic Plane. A rea- 
sonably large data set was collected and the spectrum of this 
source was measured up to energies of 1.6 TeV The differen- 
tial energy spectrum can be fitted with a power law of slope 
r = -3.1 ±0.3. The coincidence of the VHE 7-ray source 
with SNR IC 443 poses this SNR as a natural counterpart, and 
although the mechanism responsible for the high energy radi- 
ation remains yet to be clarified, a massive molecular cloud 
and OH maser emissions are located at the same sky position 
as that of MAGIC J0616+225, and suggest that a nucleonic 
origin of the VHE 7-rays is possible. 
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